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Abstract 

Approximate approach is suggested for investigation of equilibrium 
stellar models, a relativistic collapse problem and black hole formation, 
based on a Galerkin method. Some results of its simplified version - 
energetic method- are reviewed, and equation for general Gelerkin method 
are presented. 

1 Introduction 

Relativistic collapse leading to a black hole formation is a final stage of the 
evolution of the sufficiently massive stars with a mass M > (30 — 50)M©. 
To describe these collapse dynamical equations in general relativity (GR) are 
solved numerically (Baumgarte et al., f995, f996). In order to simplify the 
problem for obtaining a rapid approximate answer for such a questions, like 
collapse behaviour of different stellar masses, comparative neutrino light curves 
and their spectra, a simplified method is suggested. It is based on the classical 
Galerkin method, widely used in different mechanical and physical problems, 
connected with numerical computations (Fletcher, 1984). For the problem of 
stellar collapse this method may be considered as a generalization of the well 
known energetic method, which was applied for calculations of stellar stabil- 
ity (Zeldovich and Novikov, 1965, 1967) and initial stages of stellar collapse 
(Bisnovatyi-Kogan, 1968) in post-newtonian approximation (PN). Here we for- 
mulate the Galerkin method for the problems of relativistic collapse in GR. 
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2 Energetic method in post-newtonian approx- 
imation (PN) 

Consider first the formulation of the energetic method in PN. The main idea of 
this method is an averaging of the static or dynamic equations over the volume, 
using a given prescribed density distribution inside the star. We consider here 
only spherically symmetric stars, where the density is prescribed by polytropic 
distribution with a central density, obtained from algebraic (static) or ordinary 
differential (dynamics) equation. Consider polytropic equation of state P = 
Kp 1 . Combining the equilibrium and continuity equations, we get the equation 
for p(r) 



d 
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Where r is the Euler radial coordinate, p is the density. Transforming to di- 
mensionless quantities ( Chandrasekhar, 1939 ) 9 and £, defined as 



P 



(n + l)K i_i 



1/2 



7 



1 

n 



(2) 



where n is an polytrope index, p c is the central density, gives the equilibrium 
Lane-Emden equation 
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with the boundary conditions — 1, ^|=0 at £ = 0. The stellar boundary 
corresponds to £ = £i so that = 0. The stellar mass M expressed in terms 
of variables (0) becomes 
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and hence from (a), 
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Obviously, at n — 3, 7 = 4/3 the stellar mass is independent of p c and is 
exactly determined by the constant K in the equation of state. The stellar 
mass increases at 7 > 4/3 with increasing p c , and at 7 < 4/3 decreases. If the 
polytropic power coincides with the adiabatic power, 7 = 7 ac i, then the star is 
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stable at 7 > 4/3 and unstable at 7 < 4/3, and 7 = 4/3 corresponds to the 
boundary case and represents the indifferent equilibrium. 

Real stars are not polytropes, but the condition 7 = 4/3 is approximately 
valid at the boundary of stability if 7 is treated as an adiabatic power properly 
averaged over the star. A strict derivation of stability conditions, including GR, 
is made by use of the variational method (Harrison et al., 1965). Equating the 
first variation to zero yields the equilibrium equation, while the stability con- 
dition requires the second variation to be positive. For an isentropic polytrope 
with 7 = 4/3, p c is arbitrary, whereas the density distribution is invariant 
against homologous contraction or expansion. Let us treat these properties as 
valid also for the case where is 7 = 4/3 only in average. We then derive the equi- 
librium and stability conditions, using the simplified variational method based 
on the assumption of homology and conservation of stellar structure at den- 
sity variations (Zcldovich and Novikov, 1967), usually called energetic method. 
Write down the total energy of an instantaneously static star analogous to the 
potential energy of consevative mechanical system: 

C M C M Gmdm G 2 M 3 

e= E(p, T)dm- / 5.06 , dm = 4npr 2 dr (6) 

Jo Jo r R c 

The first term here represents the internal energy the second the new- 
tonian gravitational energy ea, and the third, £gr, a small GR correction 
{-f — ^§7? <C 1 is the small parameter) evaluated by Zeldovich and Novikov 
(1967) for the matter distribution over a n = 3 polytrope. The term contain- 
ing the Newtonian gravitational energy of equilibrium star may be explicitely 
evaluated for arbitrary polytropic equation of state (Landau and Lifshits, 1976). 

3 GM 2 

£G = -——- (?) 

For an adiabat with 7 = 7 a d, we have E — n— and for a star in equilibrium 
takes place 

n n GM 2 n-3 GM 2 

£i = -^£G^- 5-, £N = £i + £G = 7 5— > (8) 

3 5 — nit 5 — n 6 

Where £n is the total energy of a newtonian star. The total energy of a stable 
star is negative, therfore the stability requires that n < 3, 7 > 4/3. The radius 
of a polytrope is, using (||), (Q) and (||) 



R = afi 



4ttG 
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Here are used the values for a polytrope of n — 3: £1 = 6.89685, M3 = 2.01824. 



The ratio of p c to the average density p (M = ^-pR 3 ) is ^ - i:; 
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54.18. From (g), (§) and (§) we have 



CM 1 l z 

e G = -0.639 GA/ 5/3 pj /3 , g Gfl = -0.918 pg /3 . (10) 

1/3 

Only one parameter, p c or R, varies with homologous variations: 

TYl 

p = p c ip{v), v = — , <p{v) is an invariant function (11) 



and hence, the energy variations reduce to ordinary derivatives. Using (|10[)- 
( p"l| ) and taking the entropy to be constant at variations, we get from (^) the 
equilibrium condition 



° 3p 4 / 3 
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The second derivative of the energy at S = const turns into zero on the boundary 
of stability: 

d 2 £ n f M ( 4\ „ dm ln ,G 2 M 7 / 3 n 

9,-/3/ f 7 __jp__ 1 . 84 __ = o. (13) 
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We use here the thermodynamic relations 



dP 
dpi 

Equations ( |T^ ) and ( |l3| ) have been obtained by Bisnovatyi-Kogan (1966) and 
used for determining the boundary of stability for various stellar models. 



- -*«/•(§£) rf,,-*^, 7-*-(§^),. (U, 



3 Supermassive stars with a hot dark matter 

For M > 10 4 M Q the main reason of instability are GR effects. The entropy 
of such supermassive stars in critical state is so large that the pressure is de- 
termined mainly by the radiation with a small admixture of plasma, important 
for stability, but giving a very short time until the onset of instability. A com- 
mon way to overcome this instability is to consider rotating superstars, what 
may postpone the moment of collapse to 3 x 10 4 years for solid body rota- 
tion with angular momentum and mass losses (Bisnovatyi-Kogan, Zeldovich 
and Novikov, 1967), and much longer for a differentially rotating star evolv- 
ing with almost constant angular momentum (Fowler, 1966; Bisnovatyi-Kogan 
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and Ruzmaikin, 1973). Formation of supermassive stars on early stages of the 
Universe expansion, their loss of stability with subsequent collapse or explosion 
(Bisnovatyi-Kogan, 1968; Fricke, 1973; Fuller et al, 1986) could be important 
for early formation of heavy elements, observed in spectra of the most distant 
objects with red shift ~ 5, creation of perturbations for large scale structure 
formation, influence on small scale fluctuations of microwave background radia- 
tion (Peebles, 1987; Cen et al, 1993). A necessity of a presence of a dark matter 
in modern cosmological models makes it important to include it into stability 
analysis of supermassive stars. This was done by McLaughlin and Fuller (1996), 
who dealed with nonrotating superstars. The same problem for rotating super- 
stars, using energetic method, was solved by Bisnovatyi-Kogan (1998). The 
rotational effects occure to be more important for realistic choice of parameters. 



3.1 Stability analysis 

In supermassive stars with equation of state P = P r + P g = — — h pRT, 
( where TZ is a gas constant, a is a constant of a radiative energy density ) there 
is P r 3> P g due to high entropy of such stars. Besides, such stars are fully 
convective and entropy is uniform over them, so the spatial structure is well 
described by a polytropic distribution, corresponding to 7 = 4/3. The influence 
of a hot dark matter, which density does not change during perturbations, 
should be taken by account of a newtonian gravitational energy of the star in 
the dark matter potential, because GR effects of a dark matter are of a higher 
order of magnitude (McLaughlin, Fuller, 1996). For radiation dominated plasma 
there is a following expression for the adiabatic index, determining the stability 
to a collapse 

f d\o S P \ „ 4 / ft\ 4 /3 
7= l^g7j s Pa H 1 + ^J = 3 + 6' (15) 

where (3 = ^ = ^ . In the radiationaly dominated supermassive star there is a 
unique connection between its mass M and entropy per unit mass S (Zeldovich 
and Novikov, 1967) 

/_a_\3/2/3S 



M = 4 - 44 v3gJ {-a) ' < 16 > 

where a is a constant of the radiation energy density, and numerical coefficient 
is related to the polytropic density distribution with 7 = 4/3. At the point 
of a loss of stability the critical value of an average adiabatic index < 7 > in 
selfgravitating nonrotating star with account of post-newtonian corrections is 
determined by a relation (Zeldovich and Novikov, 1967) 



4 4 2 equ 4 GM 2 / 3 p 1 J 3 

<7>crs=^ + 6 GR = - + -^^-+0.99 (17) 

3 3 3 £q 3 c z 
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Here averaging is done according to (|l3|). From comparison between ( |l5| ) and 
( fl7j ) we get a well known relation for a critical central density of a supermassive 
star stabilized by plasma 

* = °- 10 gSStJ m " /2 w l8 x 1018 (^) 7/2g/cm3 ' (18) 

Here and below we consider for simplicity a pure hydrogen plasma. Newtonian 
energy of a superstar e nc i in the gravitational field of uniformly distributed dark 
matter with a density pd is written as 

l-M 

e n d = I fddm. (19) 



JO 

The gravitational potential of a uniform dark matter tfd is written as 

2tt 2 3 GM d 

where Rd is much larger then stellar radius R, and Md is a total mass of the dark 
matter halo. Stability does not depend on normalization of the gravitational 
potential so we shall omit the constant value in (fil3). It follows from ( |T9| ) and 

i 2/3 2/3 1/3 

(J) that during variations e n d ~ p c , while eqr, ~ Pc and £g ~ Pc 



(Zeldovich and Novikov, 1967). Critical value of the adiabatic index is obtained 
from (|l^), and for nonrotating superstar in presence of a dark matter the critical 
value of an average adiabatic index < 7 > crn rot is determined by 

< 7 >crnrot = < 7 >crs +Odm = 77 + 77 A f ■ 

3 3 e G |£g| 

The relation for a critical density in presence of a dark matter is obtained by 
comparison of (|l5|) and (pl|). We get 



0.99 GM y /3 = 4.1** + I = 4.1* + * fill) f ^) 3/4 . (22) 
c 2 p c 6 p c 2a V M J V3G/ v ; 

Here relations (|9|) is used for R, and J Q pr A dr = 6.95 x 10~ 4 p c i? 5 , based on 
Emdcn polytropic distribution with 7 = 4/3, was used (see Bisnovatyi-Kogan, 
1989). Relation (|§) is reduced to 

2-8 x 10- 3 ^ p 4 / 3 = 3.5 x 1(T 4 V2 p c + p d . (23) 

Solution of d23| ) is presented in Fig.l. 
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3.2 Stability of rotating stars 



Consider a rigid rotation, when its energy is a small correction to the energy of 
radiation and the energetic method is a good approach. When losses of an an- 
gular momentum during evolution are negligible we distinguish between rapidly 
rotating (RR) and slowly rotating (SR) superstars. In RR case a superstar 
reaches the state of rotational equatorial breaking before loosing its dynamical 
instability, and in SR case instability comes first. If a superstar has an angular 
momentum J, then its rotational energy e rot ~ 1.25 J 2 pc^ 3 M~ 5 / 3 , and a ra- 
tio s ro t/sGR remains constant during evolution (Bisnovatyi-Kogan, Zeldovich, 
Novikov, 1967). In presence of rotation and dark matter the critical value of 
the adiabatic index < 7 > crro t is written as 

4 2 \e GR \ - e rot s nd 
3 3 \e G \ I eel 

and the relations for determination of a critical central density, instead of (p3|), 
is written as 



2 - 8 * (£) ^ (' - S) - 35 * 10 ~ 4 (§) 1/2 * + «■ < 25 > 

As follows from (63), a superstar does dot loose its stability when e ro t > \egr\- 
This qualitative result, obtained in the post-newtonian approximation, remains 
to be valid in a strong gravitational field and reflects a presence of a limiting 
specific angular momentum au m — GM/c, so that a black holes with a Kerr 
metric may exist only at a < au m (Misner, Thorne, Wheeler, 1973). 

A RR superstar in a course of the evolution reaches instead a limit of a 
rotational instability, and equatorial mass shedding begins, leading to a loss 
of an angular momentum. Such star will loose the stability when the anuglar 
momentum will become less then the limiting value. The stage of a mass loss 
was examined by (Bisnovatyi-Kogan, Zeldovich, Novikov, 1967), where it was 
shown that this stage may last about 10 times longer, then a maximum evolution 
time to approach the rotational instability point. RR star reaches the stage of 
a rotational instability at different central densities, depending on J, but the 
ratio of rotatonal and Newtonan gravitational energy on the mass-shedding 
curve remains constant (Bisnovatyi-Kogan, Zeldovich, Novikov, 1967) [] 

e rot = 0.00725|e c |. (26) 

The energy of a rotating supertar in equilibrium in presence of a hot dark matter 
may be written as 

1 Note, that in presence of a dark matter halo mass outflow begins at rotational energy 
approximately (1 + 54p I j/p t; ) times larger, then in ( Eg ) due to additional gravity of a dark 
matter. This correction is small for a considered halo density on the mass-shedding curve at 
J > Jo (see below), and is neglected in this section. 
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£eq = -£gas + \£GIl\ - £rot + '&£nd- (27) 

In the main term for a superstar in equilibrium a relation is valid 



gas - (28) 



Taking into account (|26|), (28), we get an expression for an equilibrium energy 



along the mass-shedding curve (with variable J) 



s eq = - ( 0.00725 + f ) |e G | + \£ GR \ + 3e nd . (29) 



2, 

The curve £eq{pc) has a minimum at the central density, determined by a relation 

/ M \ 2 ^ 3 I ' M \ 
2.8 x 10- 3 f — J p 4 / 3 = 3.5 x 10- 4 f -^J p c + 5.9 x lQ- 4 p c + p d . (30) 



From comparison ( p0[ ) and (|25|) with account of (|2^) it is clear, that dynami- 
cal instability cannot occure in the minimum of the mass-shedding curve, and 
after crossing it the evolution proceeds with a substantial mass and angular mo- 
mentum losses. Central density of the superstar in the minimum of the mass- 
shedding curve ( p9| ) with and without dark matter are represented in the fig.1. 

Let us find a parameters of a superstar, at which its critical state is situated 
on the mass-shedding curve. These parameters satisfy sumultanously the rela- 



tions (25) and (|26|), what leads to the equation for determination of a central 
density in the form 



2.8 x 10- 3 (^) ^ P 4 J 3 = 3-5 x 10- 4 (Ml) V2 Pc + 12 x lQ- 4 p c + Pd . (31) 



The relation ( |26| ) is used for determination of an angular momentum of the 
superstar J = Jq with p c from (30) and J = J\ < Jo with p c from (|3l|). Solution 



of this equation is also given in fig.l. As may be seen from fig.l, the stabilizing 
effect of rotation on the mass-shedding curve at J = J\ is more important, then 
stabilization by a hot dark matter, which influence decreases at increasing of 
a central density. With account of a longer state of the evolution with mass 
loss until reaching the point of the loss of stability at larger J, this conclusion 
becomes even stronger. 



4 Collapse and explosions of supermassive stars 

To study dynamical processes by the energetic method we use, instead of the 
energy variation, the energy conservation law in the form 
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Figure 1: The correction terms Sqr and 5 dm + $GRi the quantities /3/6 (line c), 
l3/6+{s ro t/£N)sh/3 (line b), and 2(3/6+ (s ro t/sN)sh/^ (line a), as functions of 
the central density of a supermassive star with M = 1O 6 M0, and dark matter 
density of 10 -5 g/cm 3 . The instability points for nonrotating star occure at 
intersection of the correction term curves with the line c. Mass shedding in 
the stable star with angular momentum Jo (see text) occures at intersection of 
correction term curves with the line b, and critical point on the mass-shedding 
curve is determined by a corresponding intersection with the line a. 
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d(e in / dQdm = dS Tdm. (32) 

Jo Jo 

Here Sk — \ f v 2 dm is the kinetic energy of the superstar. Using thermody- 
namic relation we get 



r M r RI f QE\ f M f dE\ 

dei n — d E(p, S)dm = / ( — — I dp dm + / ( —— I dS dm 
Jo Jo \dpJs Jo \dS J P 



— dp c I — dm + dS I Tdm. (33) 



la P 2 dp c 

From the mass conservation law, in presence of homologycal motion with the 
fixed density distribution in space (pd|), we obtain a space velocity distribution. 
Considering r as a Lagrangian radius, corresponding to the mass m < M, 
r(M) = R, we may write 

m = A{£)p c r 3 , M = A oPc R 3 , = A Q = -g-M n , 

dRc +3 Pcdr = dR 1+3 PcdR = 

dt r dt ' dt R dt K ' 

Using definitions of the velocity in Lagrangian coordinates, we get its homolog- 
ical space distribution 

dr dR r 

v = Tr = ~dt' v = Vr r- (35) 

With account of the equality § Q l {r 2 / R 2 )ip(v)dv = 0.108, we get a dynamical 
equation in the form 

0.597 {P :\ > + 0.639 Gp 2 J 3 + 1.84 - , p c (36) 
dt z <r 

-3M- 2 / 3 p: 2 / 3 C P^- =0. 
Jo 

On dynamical stages of the superstar evolution, after its loss of stability the 
radiative transfer losses (— Q r ) are less important, then neutrino losses (— Q v ) 
and energy production in nuclear reactions Q n . The equation determining en- 
tropy changes is averaged over the star with a weight function (T 4 ~ Ep) for 
radiation dominated superstar, which entropy is taken homogenous 



dl f 1 

— / pT^dv = < T 4 Q n > - < T 4 Q v > - < T 4 Q r > . 
dt J 



(37) 
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Nuclear reactions of pp and CNO hydrogen burning and 3a helium burning have 
been included in the calculations of Bisnovatyi-Kogan (1968). The equaions 
(|36"|) , d37|) , together with equations for averaged composition of hydrogen X, 
helium Y, and equation of state with account of e + e~ pairs creation and rela- 
tivists relations for electrons, had been solved numerically. Primodial chemical 
composition with only hydrogen and helium was taken as initial condition. In 
the process of contraction after a loss of stability, 3a reaction produces 12 C, 
which initiates a CNO hydrogen burning, pp reaction remaining always unim- 
portant. It was obtained that in stars with M < 1.5 x 10 5 M Q collapse is 
reversed, and they explode, enriching the intergalactic and interstellar gas with 
heavy elements. Such explosions could happen on stages, preceding the epoch of 
a galaxy formation. That is one of the way to explain high metallicity in the dis- 
tant quasars and intergalactic gas in galaxy clusters. Similar calculations made 
for rotating superstars, and for normal (solar) composition shift the boundary 
between collapsing and exploding superstars to higher masses (Fricke,1973). 

5 Formulation of Galerkin method 

The Galerkin method allows to find an approximate solution of differential equa- 
tions. In this method the solution of the partial or ordinary differential equation 
is represented in the form (Fletcher, 1984) 

N 

u = up +}^ajipi, (38) 
i=l 

where ifi are known analytical functions, and coefficients (or functions) a% are to 
be determined. Functions (pi are assumed to be not necessary orthogonal. When 
substituted to the original differential equation the coefficients a, will satisfy a 
set of algebraic or ordinary differential equations, obtained by minimizing the 
corresponding functional. 

Consider how the Galerkin method can be applied to the problem of the 
relativistic star collapse. The functional which should be minimized, is a total 
stellar energy. It contains a kinetic energy ( |32| ) , which in PN is written explicitly. 
In GR the total energy cannot be splittcd into different components, and the 
collapsing body should be described by means of a stress energy tensor X^. 
Equating to zero the variation of the total energy gives dynamical equations of 
the collapse. In the process of contraction the collapse accelerates and finally 
approaches a free fall which is described by a Tolman solution (Landau and 
Lifshits, 1962). The approximate method should be formulated in such a way, 
that its solution approaches the Tolman one at final stages. 

We may use the Galerkin method to approximate the system of equations in 
partial derivatives by the system of ordinary differential equations. The density, 
as a function of a lagrangian coordinate a and time t, p(a, t) is represented as a 
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sum of ati(t)(pi(a). Since tp(a) is invariant under the lagrangian transformations, 
the partial time derivative of this product will transform to ^-ipi{a). Thus, the 
Galerkin method allows to separate variables t, and a. The same separation has 
to be done for a velocity function v(a, t). 

5.1 Stability conditions in PN 

In PN approximation the density in Galerkin method is written as a following 
sum 



N N 



au(t)<pi(a) , where p c = ^ oti(t)cp(Q) (39) 



As for a function ipo(a), it is convenient to take corresponding Emden profile for 
one of polytropic indices. For other functions we may chose tp^ = cos K Tta, 
which have increasing number of nodes. 

Then satisfaction of the boundary conditions: 

<Pi(A) = 0, A = a(R); p,(0) = 1 (40) 

will be provided. The minimization of the energy functional (^|) for finding an 
equilibrium model is reduced to zero partial derivatives 

leading in the static case of constant cti to a set of N algebraic equations for 
finding equilibrium ctl . Stability of a model is found from an evaluation of the 
second variation <5 2 e, which in the energetic method is reduced to the algebraic 
equation ( |l3| ) . In the Galerkin method with several scaling functions tpi (a) , the 
second variation 8 2 e is represented by a quadratic form 

N f) 2 

S 2 e = J2^- Sa i Sa k, (42) 

l,K 

In the energetic method the condition of stability reduces to one equation 
d 2 e/da 2 > 0, given in (p^). In the Galerkin method the stability is related 
to positive definiteness of the quadratic form (pl2|), what is provided (Smirnov, 
1958) by the positiveness of the determinant 

-II > 0, (43) 



dondcik 

and all its main minors. Remind, that main minors are determinants, obtaind 
from the main determinant ( fi"3| ) after eliminating lines and columns, intersecting 
on the main diagonal at i = k. For two functions in (EM) the positiveness 
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of the main determinant (|43|), and two partial derivatives d 2 e/da\ > and 
d 2 e/da\ > are enough for stellar stability. Loss of stability happens close 
before the point where the main determinant, or one of its main minors becomes 
zero. 

In approximate presentation of the trial function in the Galerkin method, 
the minimal value of the second energy variation is larger, then its value for a 
real trial function. So zero values of the determinant ([43|), or one of its main 
minors, guarantees the onset of instability. Their positiveness is not an exact 
guarantee of the stability, but comparison of the energetic method with an exact 
stability analysis shows a good presicion of this approximate approach in most 
realistic cases (Bisnovatyi-Kogan, 1989). Energetic method corresponds to a 
homologeous trial function for displacement Sr ~ r. In the Galerkin method the 
trial function may be determined with a better precision. In fact, the coefficients 
8on for the trial function of a density 

N 

5p = y ^2 / 8a i ip i {a) (44) 

i 

are determined as an eigenvector of a set of uniform linear equations 

d 2 - 



Sa k = XpSoti. (45) 

The eigenvector Saf is used for obtaining an approximate eigenfunction (Q), 
and eigenvalues A p are related to the square eigenfrequcncies of the stellar model. 
The positive definiteness of the quadratic form ([l|]) coincides with the positive- 
ness of all eigenvalues X p . So, Galerkin method permits to investigate a stability 
by finding approximate eigenvalues and eigenfunctions of a linear set of algebraic 
equations, instead of finding the same values from a second order differential 
equation exactly. 

5.2 Relations for total energy and barion number in GR 



In GR a barion number density n is used in ( p9[ ) instead of the mass density 
p, which in this case is not presenting itself a conserved value. To obtain GR 
expression for the total energy of a spherically symmetric radially moving body 
take a stress-energy tensor for the perfect gas 

Tij = (P + ne)uiUj - Pg%j, (46) 

where e is an internal energy per baryon, Ui is a four velocity of the matter. 
Inside the spherically symmetric star the Schwarzschild type metric is used 

gij = e^(?dt 2 - e X dr 2 -r 2 dVl 2 , (47) 
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where dQ 2 = dO 2 + sin 2 Odip 2 , and A is expressed as a functional of Tq (Landau 
and Lifshits, 1962) 

r 

e > = (l _ e - 4tt / T V dr. (48) 

c 4 r J 


To avoid confusion, we determine r and i to be independent Euler variables 
and v to be a velocity field of the medium. A number of baryons a inside a 
Lagrangian radius r(a) is considered as a Lagrangian independent coordinate 

r nr 2 

a = An ; dr. (49) 

_ 2Ge 

Then, the world line of each shell is r(i, a). There exist a set of transformations 
of the Euler type metric (fl7|) to its Lagrangian representation: 

9tj = e 2v c 2 dT 2 - e 2A da 2 - r 2 dn 2 , (50) 

where r scales time of the co-moving observer. We determine a physical velocity 
of the shell as a displacement dl, measured by the Euler rest observer, with 
respect to Euler time, measured by the same observer: 



0/2 f g r 



c 4 r 



(51) 



Where (j> to be determined from ( f47| ) . In a case of a test particle it gives a well 
known result for the free falling particle in a Schwarzschild metric 

2GM. 1 dr 2GM x . . dr r 

v = [l — ) -jr = (l — ) ^ r = 17 =w - ( 52 ) 

c r at c r ar 

Covariant components v a and are given by the following relations (Landau 
and Lifshits, 1962) 

v a = Japv 13 , v 2 = v a v a . (53) 

In a case of a diagonal metric, the dimensionless 3 - metric tensor is connected 
with the 4 - metric tensor as — —g a p- The mass-energy functional for the 
spherically symmetrical relativistic star is written as in (fl§| ) 

R 

E = e{R) = 4tt J T °r 2 dr. (54) 
o 

Four velocity of the fluid can be expressed in terms of the tree velocity (Landau 
and Lifshits, 1962) 
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Then, T® component of the stress-energy tensor ( |46| ) is written as 

T ° = • (56) 

From (E^,p9) making use of (p5^) we get 



de ne + P^ I 2Ge 



On 1-4 V 1 c 4 r' (57) 



— - 1 /i 2Ge 

5.3 Velocity distribution over the star in GR 

Assuming the profile law for n, let us obtain the profile law for the velocity. To 
separate variables a and t we need a relation between the pertubations of the 
boundary (dr/dt) a=N , and velocity of the inner shell (dr/dt) a . Differentiating 
©gives 



dr _ 1 -2Ge/c 4 r + 2Ger/c 4 r 2 V 1- ^/" 



da 87mr 2 2Ge 27mr 2 \ 2n 

c 4 r 



(59) 



what contains terms with e. After differentiation, equation ( |57| ) contains a 
second order derivative r. Going on this iterative process, we obtain an infinite 
system of differential equations: 



da (1-f) + c 2 (l-^) 2 

1 (nc + P4) (-2Ge/c 4 r + 2Gef/c 4 r 2 ) 



2 (i-SWi-^) 



<3e 
9a 



(60) 

(61) 
(62) 



completed by a similar set of equations obtained from differentiation of (p9[). 
In PN the velocity distribution (p5[) is found from the continuity equation, and 
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than is used in derivation of the equation of motion (pq) . In full GR description 



the infinite set of equations (60)-(p2[) should be used for obtaining the velocity 
distribution, instead of (|35|). Such complication is connected with an influence 
of the 'kinetic' energy on the geometry of the star interior. After substituting 
of expansion on n, similar to (|3^), we obtain ordinary differential equations 
instead of the partial ones. In practice the system ( |60| ) - (|62| ) has to be cutted 
off on some step. Final stages of the gravitational collapse are close to a free 
fall of a dust spherical cloud. Thus using e from the Tolman solution should be 
good for late stages of the collapse. For the zero-pressure dust collapse without 
thermal processes the energy inside a given lagrangian shell does not change, so 



in fl59| ) the relation e(a) = (de/dt) a — could be approximately used, making 
unnesessary to solve the equations ( |60|) - ( |62] ) for obtaining the profile v(a). 

5.4 Equations of motion 

The equation determining the change of the entropy of the star is written similar 
to (|37|). In that case the entropy distribution over the star should not change 
its form in the process of the collapse, or be uniform. The integration should 
be done over the lagrangian coordinate da — 4 T" r ,^. r e . To derive equation 

of motion we should, like in (J3q), to take derivatives from all thermodynamic 
functions at constant entropy distribution, as follows from (|32|) , (|33|) . 

When only one term is is left in (|39|), as in the energetic method, the equation 
( |60|) is used also for derivation of the equation of motion. In the case of several 
terms in ([}9|), the equation ( |60|) should be written for each coefficient at, taking 
constant other on, i 7^ k. This the number of equations will be equal to the 
number of unknown functions a, . This procedure is equivalent to a minimization 
of the energy functional over all coefficients oii(t), or over all values of on for the 
equilibrium solution. 

It is important to stress, that when using (^0|) - ([32]) for derivation of the 
velocity profile, all coefficients a,(£) must be differentiated together. 



6 Conclusion 

We suggest an approximate method for solving the problem of stellar stability, 
and dynamical stages of relativistic stellar collapse, using Galerkin method. It 
should be mentioned, that this paper, together with the general review of the 
contemporary state of art of the problem contains the description of the essen- 
tially new results. The stability analysis of the rotating stars with the presence 
of the surrounding hot dark matter background, as well as the application of 
the energetic method together with the Galerkin method to the problem of 
relativistic gravitational collapse are such new results. 

In the case of stability investigation we get a set of algebraic equations and 
search of the eigenvalues of matrics, instead of solution of eigenvalue problem 
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for differential equation of a second order. In the case of collapse set of ordinary 
differential equations should be solved, instead of partial ones. It is evident, that 
solution of ordinary equations is much simpler, what could permit to investigate 
a wide set of equation of states, stellar masses, and to get different neutrino light 
curves during a black hole formation. 

Several examples of application of the simplest version of Galcrkin method, 
a well known energetic method, are presented. A simple solution obtained by 
this method for complicated problems show explicitely the advantage of such 
approximate approach. 
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